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Abstract 

We report about the energy and momentum resolved optical response of black phosphorus (BP) in 
its bulk form. Along the armchair direction of the puckered layers we find a highly dispersive mode 
that is strongly suppressed in the perpendicular (zig-zag) direction. This mode emerges out of the 
single-particle continuum for finite values of momentum and is therefore interpreted as an exciton. 
We argue that this exciton, which has already been predicted theoretically for phosphorene - the 
monolayer form of BP - can be detected by conventional optical spectroscopy in the two-dimensional 
case and might pave the way for optoelectronic applications of this emerging material. 
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Among the different materials which are currently under investigation in the “post¬ 
graphene” era, black phosphorus (BP) in its monolayer form phosphorene is of particular 
interest [1]. Due to its intrinsic band-gap of about 0.3 eV in the bulk and much higher values 
for the monolayer [2, 3], devices built from BP show substantially increased performance 
compared to pristine graphene in terms of e.g. on/off ratios [3, 4], 

In addition, very high mobility values (up to 10 4 cm 2 V~ 4 s _1 ) have been predicted 
theoretically for phosphorene [5] and at least to some extend confirmed experimentally [4], 
What is particularly remarkable is a strong in-plane anisotropy of the transport and optical 
properties [5-9] that is, however, well known already for the bulk [10, 11]. Even more, it has 
been predicted that phosphorene allows for the fabrication of fullerene- and nanotube-like 
aggregates [12] with highly tunable electronic properties. Together with the known tendency 
to undergo several structural phase transitions under pressure - which eventually results in a 
superconducting ground state with T c & 10 K [13] - this promotes BP to a versatile platform 
for possible applications. 

While a thorough understanding of the charge response is of utmost importance for all 
optical applications, there is a surprising lack of corresponding data, in particular in the 
near-infrared to visible regime. 

Here, by employing electron energy-loss spectroscopy (EELS) in transmission, we bridge 
this gap. We find strongly anisotropic optical properties which are the result of an excitonic 
state that is strongly polarized and - at least in the bulk of BP - exists only away from the 
center of the Brillouin zone (BZ). 

EELS in transmission is a bulk-sensitive scattering technique whose cross-section is propor¬ 
tional to the so called loss-function L(q } co) = Im(—l/e(g,cu)), with e(q,uj) the momentum- 
and energy-resolved dielectric function [14]. It has been successfully applied to investigate 
collective charge modes in numerous condensed-matter systems (see [15, 16] for an overview). 
In particular, it has proven useful for the investigation of excitonic states in organic semi¬ 
conductors [17-19]. In addition it is possible to retrieve the complete dielectric function 
e(q,uj) = 6i(q,co) + ie 2 (q,uj) via a Kramers-Kronig analysis (KKA) of the measured data. 

For the present experiments thin films have been prepared via exfoliation from a single 
crystal purchased from “2d Semiconductors Inc.” [20]. The films (d ~ 100 — 200 nm) were 
put on standard electon-microscopy grids and then transferred to the spectrometer where 
they have been aligned in-situ with electron diffraction showing the high quality of our 
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samples and allowing for polarization dependent investigations along well defined directions 
within the puckered layers of BP. The measurements were carried out using a purpose-built 
transmission electron energy-loss spectrometer [21] with a primary electron energy of 172 keV 
and energy and momentum resolutions of A E = 80 meV and A q = 0.035 A respectively 
and at temperatures of T « 20 K to minimize thermal broadening. All discussed features 
remain qualitatively equivalent at room temperature. 

First-principles calculations have been performed with the FPLO package [22], We 
employed the crystal structure given in [23] except for the in-plane lattice constants which 
we determined experimentally from electron diffraction to be 4. 363A along the armchair and 
3.206A along the zig-zag direction, respectively. The total density was converged on a grid 
of 12 x 12 x 12 irreducible /e-points. The generalized gradient approximation (GGA) for 
the exchange correlation potential has been used as parametrized by Perdew, Burke and 
Ernzerhof [24], 

Under ambient conditions BP in its bulk form adopts an orthorhombic lattice with 
puckered layers spanned by the x- and y -axes along which the atoms are arranged in armchair 
and zig-zag coordination, respectively. This symmetry allows for two independent in-plane 
components of the dielectric tensor e xy (q,u). 

From Fig. 1 it is clear that this translates directly to a strongly anisotropic EELS response. 
Along the TX (armchair) direction we find a low-energy mode that resides at about E = 0.6 eV, 
slightly above the gap edge of Eg = 0.3 eV [25]. When changing the polarization direction 
from the TX- to the TY (zig-zag)-direction within the puckered layers, this mode gradually 
shifts to lower energy and loses strength as can be seen from Fig. 1(b). 

Considering a simple Lorentz picture for the observed structure it is important to note 
that a peak in the EELS response occurs not at the energy uq of the electronic transition but 
at uj = y/cjg + Wp where ui p is a measure of the spectral weight acquired by this excitation 
[26]. It can be shown that in the limit of —>■ 1, where is the background dielectric 
constant, L(Cj) ~ u p /uq which means that a concomitant lowering of intensity and energetic 
position of a transition observed in EELS can only be reconciled by a diminished spectral 
weight Up. Consequently, we attribute the behavior of the measured EELS intensity shown in 
Fig. 1 to an angular dependent reduction of spectral weight when changing the polarization 
from TX to TY. 

Moving away from the center of the BZ this low-energy mode shows a strong momentum 
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Figure 1: (Color online) (a) EELS intensity along the two high-symmetry directions within 
the puckered layers of BP. (b) Low-energy polarization-map showing the smooth angular 
dependence of the EELS intensity within the puckered layers of BP. (c) Calculated 
single-particle band structure along TX. Data in (a) and (b) have been measured with a 

fixed momentum transfer of q — 0.1 A 


dependence (see Fig. 2(a)) which hints at a highly dispersive and therefore delocalized 
excitation. This is further corroborated by Fig. 2(b) where we show the dispersion of this 
feature as evaluated by tracking the steepest slope on the low-energy side. We emphasize that 
along the TY direction there is no dispersion and the spectrum remains largely unaffected 
in the corresponding energy window when increasing the momentum transfer. In addition, 
Fig. 2(b) contains the size of the single-particle gap EQ FT (q ) = Ecb(q) — Evb{ T) as derived 
from our DFT calculations. 

As can be seen, the feature under investigations shows a bandwidth of more than 1.5 eV 
within the accessible momentum window. More importantly, when fitted to a quadratic 
function, its onset Eq = E(q = 0) = 0.473 eV lies in close proximity to the single-particle 
gap. Note that we cannot measure at lower g-values because of the increasing contribution 
of surface scattering for q —> 0 [21], Comparing the bandgap extracted from our DFT 
calculations to the measured EELS data, there is clear evidence that the dispersion of 
the low-energy mode seen in EELS lies below the single-particle gap at finite q which is 
suggestive of a bound, i.e. excitonic state. The fit to our data shown in Fig. 2(b) yields 
an effective exciton mass of M* = h 2 /2A ~ 1.6 m e with the free-electron mass m e . This is 
in reasonable agreement with the value of = rn* h (T) + ■ m* e (T ) ~ 1.1 m e obtained 
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Figure 2: (Color online) (a) The momentum evolution of the measured EELS intensity for 
momentum transfers along the armchair direction, (b) the dispersion of the low-energy 
feature as extracted from the low-energy edges of the measured intensities. The solid line 
corresponds to a quadratic fit with the values: A = 2.412 eV A 2 and E 0 = 0.473 eV and the 
dashed line to the single-particle gap derived from the band-structure calculation shown in 
Fig. 1. (c) The evolution of the exciton binding energy Eff as a function of momentum 

transfer along the armchair direction. 


within the effective-mass approximation for the exciton from our DFT calculations. As 
shown in Fig. 2(c), the binding energy defined as E |f = EQ FT {q) — E(q) of this exciton has 
a well pronounced maximum of about 0.1 eV for g max ~ 0.23 7r/c. We emphasize that as 
soon as E e g < 0 the exciton ceases to exist as it will become over-damped by single-particle 
excitations. As the exciton binding energy E e g(q) ~ 1 /e(q) and e(q) strongly decreases for 
finite q [27-29] the increase of E^(q) for q < q max seen in Fig. 2(c) is a result of the less 
effective screening of the Coulomb interaction for smaller wavelengths. This effect is, however, 
overcompensated by the merging of the exciton dispersion with the single particle continuum 
which results in the decay of the exciton into unbound electron-hole pairs, thereby producing 
an increasing broadening of the spectral features for higher momentum values. 

In order to obtain the full optical response, the measured EELS data have been corrected 
for contributions from the quasi-elastic line and multiple-scattering. Subsequently, the signal 
was subject to a KKA and normalized to yield €\{uj = 0) according to [30]. 

The resulting curves for the real and imaginary parts of e{uS) are shown in Fig. 3. As 
expected already from the plain EELS data shown in Fig. 1(a), along the armchair direction 
there is a well pronounced mode right above the bandgap which is invisible for polarizations 
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Figure 3: (Color online) 6\ (a) and e 2 (b) along the two high-symmetry directions within the 
puckered layers of BP as derived from a KKA of the corresponding EELS data. The insets 

show the low-energy region. 

along the zig-zag direction. The high values of seen for small momentum transfers in 
Fig. 3(a), which inevitably have to decrease for q > 0 [27-29], support our interpretation of 
an increasing E e £ (see Fig. 2(c)) in terms of weaker screening at least for q < q max . 

Strongly anisotropic optical properties have also been predicted for phosphorene by ab- 
initio calculations [2, 5]. In particular, it has been shown that phosphorene hosts an excitonic 
state which is strongly delocalized along the armchair direction but confined to a small region 
along the zig-zag direction [2, 31]. This is consistent with our interpretation given above 
for the polarization dependence of the EELS signal in the bulk. The large extension of the 
exciton along the armchair direction produces a large dipole moment and hence spectral 
weight along this direction, yielding the peak seen in Fig. 1(a) and (b). Conversely, changing 
the polarization towards FY will progressively project out a smaller dipole moment that 
lowers the intensity (and in EELS also energy) due to the small extension of the exciton 
along the zig-zag direction. 

Lowering the dimensionality from the bulk to the two-dimensional phosphorene will cause 
smaller interlayer screening resulting in a larger value of the single-particle gap [2, 5] as well 
as an enhanced excitonic binding energy due to the dependence of the dielectric function 
on the gap size [27, 28]. Recent theoretical approaches predict excitonic binding energies 
of about 0.8 eV in monolayer phosphorene [2, 31]. Therefore, in phosphorene, the exciton 
might eventually be pushed out of the single-particle continuum already at the T-point 
and produce strong optical absorption along the armchair direction. This emergence of the 
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exciton out of the continuum should then be observable as a function of layer-thickness by 
conventional optical means. Indeed, recent photoluminescence data strongly argue in favor 
of layer-dependent excitonic absorption [9]. 

From investigations on organic semiconductors it is well-known that delocalized, i.e. 
strongly dispersing excitonic modes contribute significantly to the mobility-values via the 
exciton diffusion coefficient [32, 33]. Our data therefore may provide a route to reconcile the 
high mobility values and their anisotropy predicted and reported in transport measurements 
on phosphorene [4-9] with the microscopic fate of the excitonic states in this monolayer 
material. 

In summary, we investigated the energy- and momentum-dependent optical response in 
the bulk of black phosphorus. We fold strong anisotropies within the puckered layers with 
the armchair direction showing much stronger absorption in the infrared and visible part of 
the spectrum. In addition we identified an excitonic mode for finite values of momentum 
transfer with an effective mass of about 1.6m e . This mode can be probed in phosphorene 
by conventional optical spectroscopy and might be interesting for future optoelectronic 
applications. 

We highly appreciate technical support from R. Hiibel, S. Leger and M. Naumann. 
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